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1 Introduction 

A basic strategy to obtain hypersurfaces with prescribed mean curvature 
in Euclidean space is to describe them non-parametrically as solutions of a 
Dirichlet problem for a certain quasilinear elliptic PDE. The solutions are 
then graphs over domains in totally geodesic hypersurfaces of the ambient 
space. Classical references on the subject are 0, and [16J. Existence 
results were obtained also for curved space forms in several formulations 
as one may consult 0, 0/ (H, O, EH and d. Recently the cases 



of Riemannian and warped products also deserved research efforts. For 
instance, see Q, @ and UTTI . 

The mere possibility of applying similar analytical and geometrical 
tools in all contexts seen above indicates that one must search for gen- 
eral existence results in a large class of Riemannian spaces. The presence 
of an isometric or conformal Killing vector field playing a major role is 
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certainly the distinguishing feature for all the particular ambient spaces 
aforementioned. 

In this paper, we deal with Riemannian spaces endowed with a Rilling 
vector field and define a notion of Rilling graph in these general ambi- 
ents. We are then able to solve the corresponding Dirichlet problem for 
prescribed mean curvature under hypothesis involving domain data and 
the Ricci curvature. Since the ambient metric is indeed warped in a sense 
we made precise later, the present article should be considered as an exten- 
sion and generalization of results proved in [6J. In that sense, it is worth to 
mention that key arguments in flfj do not longer work in the more general 
framework considered here. 

The Ricci curvature naturally arises in all apriori estimates we made 
since they are based on comparison of geometric data. By its turn, ge- 
ometric comparison results follow from the classical Jacobi and Ricatti 
equations. The last one is used in Section 3 for controlling the extrinsic 
geometry of the barriers. 

We now explain more precisely the framework we are considering. Let 
M be a (n + 1)— dimensional Riemannian manifold endowed with a non- 
singular Rilling vector field Y. We assume that the distribution orthogonal 
to Y is integrable. Then, the leaves are easily seen to be totally geodesic 
hypersurfaces. Let P be a fixed integral leaf and assume that the flow lines 
of the flux $:IxP->M generated by Y are complete. Given a bounded 
domain Q in P, the Rilling graph E associated to a function u on Vt is the 
hypersurface 

E = {q = *(it(p),p) :peCl}. 

Our results assure the existence of Rilling graphs with prescribed mean 
curvature H and boundary data (p. Here, the functions H and 4> are defined 
respectively on and Y where Y = OVl. The problem of existence of such 
Rilling graphs is formulated in terms of a Dirichlet problem for a diver- 
gence form elliptic PDE. The barriers we used for estimating height and 
gradient of u are the Rilling cylinders. The Rilling cylinder K over Y is 
ruled by the flow lines of Y through Y. Therefore, we have 

K = {q = #(s,p) : p e T}. 

In the sequel, the mean curvature of K pointing inward is denoted by H cy \. 

We are now in condition to state the theorems proved in this paper. We 
refer to Section 2 for the convention we used for the Ricci tensor. 
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Theorem 1. Let VI c F be a bounded domain with C 2 ' a boundary F. Sup-pose 
that H cy i > and that 

inf RicM > — n inf H 2 A . 

m r y 

Let H e C a {yt) and <fi e C 2 ' Q (T) be given. If sup Q < inf r i? C yi, f^en ffoere 
exfsfs fl unique function u e C 2 ' a (fi) satisfying u\r — <t> whose Killing graph has 
mean curvature H. 

If the Ricci tensor does not satisfy the assumption given in the above 
result, then we may use comparison theorems with geodesic spheres as 
barriers in order to prove the following result. In this situation, we must 
impose certain condition either on the function H or in the size of the do- 
main Q. 

Theorem 2. Let Q c P be a bounded domain with C 2,a boundary F contained 
in a normal geodesic disk with radius r . Suppose that inf M Ric M > —(n-l)kfor 
some positive constant k. Let H e C a {Vt) and <p e C 2,a (<9f2) be given. Assume 
that H cy i > 0, ifozi sup n \H\ < inf r # C yi #w*f 

1 x sup n |#| 

r < — p= coth 



V^ Vk 

Then there exists a unique function u e C 2 ' a (f2) satisfying u\r = </) whose Killing 
graph has mean curvature H. 

If the metric induced on P is rotationally symmetric then we may take 
certain constant mean curvature spheres as barriers. These spheres are ro- 
tationally invariant hypersurfaces and their qualitative aspect is described 
by a flux formula. In what follows, we denote r = dist(p , • ) for some 
point p E P. 

Theorem 3. Suppose that the induced metric in P is of the form 

dr 2 + Z\r)d6 2 

for a given function £, where d6 2 is the canonical metric on the unit sphere §™ _1 . 
Let Q c P be a bounded domain with C 2,a boundary F contained in a normal 
geodesic disk with radius r . Let H e C a (£l) and e C 2,a (dVt) be given. 
Assume that H cy \ > 0, that sup n \H\ < inf r H cy \ and that 

sup |/i I < 



n 1 ^VWe-H^dr 
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where g = \Y\ 2 . Then there exists a unique function u e C 2,a (Vl) satisfying 
u\t = 4> whose Killing graph has mean curvature H. 

We emphasize that in the statements and proofs of the theorems we 
may replace M (after passing to the universal cover if necessary) by the 
solid cylinder \J/(M x Vt) whose boundary is K. Moreover, the hypothe- 
sis in Theorems 2 and 3 may be understood either as restrictions on the 
size of the domains for arbitrary mean curvature functions, or as upper 
bounds on \H\ in the case of arbitrarily large domains. Finally, we remark 
that standard regularity theorems imply that the results remain true for 
continuous boundary data. 

If the assumption sup c \H\ < inf r H CJ \ fails at some point, we do not 
show here that the our results are no longer true in the sense that there is 
a boundary data <p for which no solution exists. Nevertheless, this is well 
known in some cases, namely, the Euclidean space for standard graphs 
and has also been proved for two types of Killing graphs in the hyperbolic 
space in [8 J and llTil . 

This paper is organized as follows. In Section 2, we present the basic 
geometric structure of the ambient spaces we are dealing with, including 
calculations concerning the mean curvature of the Killing cylinders and of 
the Killing graphs. In Section 3, we deduce the height estimates. Sections 
4 and 5 are devoted, respectively, to boundary and interior gradient esti- 
mates. The last section presents the proof of the theorems following the 
well-known continuity method. An appendix contains the sketched proof 
of the flux formula. 

We point out that ambient spaces with a Killing vector field correspond 
in Lorentzian setting to the important notion of stationary space-times on 
which the metric tensor is time-independent. Prescribed mean curvature 
hypersurfaces work in this context as Cauchy hypersurfaces for the initial 
value formulation for the Cauchy problem in General Relativity. Thus, it 
seems interesting to investigate the Lorentzian analogues to our existence 
results. 

The case of non-integrable orthogonal distributions corresponds to Rie- 
mannian submersions other than the simple ones associated to warped 
products. The particular situation of the three-dimensional Heisenberg Lie 
group with a left-invariant metric was treated in Q. Two of the authors of 
the pressent paper address in [5 j the question for Riemannian submersions 
with uni-dimensional totally geodesic fibers. The case of conformal vector 
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fields was studied by one of the authors in joint work with E Andrade in 
0. 

2 Killing graphs 

Let M be a connected (n+ 1) —dimensional Riemannian manifold endowed 
with a non-singular Killing vector field Y. We denote the metric and the 
Riemannian connection in M by (• , •) and V, respectively We assume that 
the flow lines of Y are complete and that the distribution 

pe M ^ {v G T P M : (v, Y) = 0} 

is integrable. Then, it is easy to verify that the integral leaves are totally 
geodesic hypersurfaces. 

Let P be such an integral leaf. The flux : R x P — > M generated by Y 
takes isometrically P to the leaves P s = ^(P), where ty s = ^(s, ■ ). Given 
local coordinates x 1 , . . . ,x n for P, then s, x 1 , . . . , x n are local coordinates for 
M defined by 

g G M (s, x 1 , . . . , x n ) if q = ^(s,p), 

where p G P is the point with coordinates x 1 , . . . , x n . The corresponding 
coordinate vector fields are 

d a = ^-*{s,p)=Y(*(s,p)) 
ds 

and 

d 

di(q) = g^^^P) = *sMdi(p)- 
The ambient metric in terms of these coordinates has components 

5-00 = (d a , d s ) = g 2 , g 0i = (d s , 9<) = 

and 

9ij = (^sA^sA) = AAj) = cry, 

where cr^ are the components of the metric in P in terms of the coordinates 
x l . Observe that the components of the metric do not depend on s. The 
gradient of the function s is 

Vs = g 00 d s = \Y\- 2 Y =: fY. 
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Since the flow lines of Y have constant geodesic curvature and parallel 
curvature vector it is a standard fact that the solid cylinder f (1 x 0) has a 
warped product structure as 

whose metric is 

<Jij dxMx- 7 + g 2 ds 2 . 

In fact, this is the setting considered in 10. Notice that the relation between 
the Ricci curvatures of M and P is determined by g = \Y\. 

Given a bounded C 2,a domain Vt on P and a smooth function u on VL, 
we define the associated Killing graph £ by 

£ = {q = *(u(p),p) : p £ f2}. 

We may think of E as the locus 

$(s,p) = s — it(p) = 0, 

where u(p) = u(s,p). An orientation for £ is given at q = ty(u(p),p) by 

= f d s -y u(p) yu{p), 

where V denotes the Riemannian connection in P and 

Vu = a %:) Uidj = u J dj 
is the gradient relatively to P. Then 

iV = i ? V$ = i ? (/S.-*.V«) (1) 

defines a unit normal vector field along £, where 

jy 2 = /+ |Vu| 2 
and ^ is a shorthand notation for *& u (p)*- 
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2.1 Killing cylinder 



The Killing cylinder over T = dQis the surface ruled by the flow lines of Y 
given by 

K = {V(s,p) : s e R, p e r}. 

If s 1 , . . . , s™^ 1 are local coordinates for T, then s, s 1 , . . . , s™ _1 are local coor- 
dinates for K. Let 77 be the unit inward normal vector along T as a sub- 
manifold of P. We equally denote by rj the unit normal vector field ^ S ,V 
along K. Thus, we have 

(r],d s ) = = (77, di). 
Since 77 and <9j are tangent to the totally geodesic leaves P s/ it results that 

Hence d s is a principal direction of K, and the corresponding principal 
curvature is the geodesic curvature 

« = f(^a s d s ,v) 

of the flow lines through T. 

In the sequel, we deduce some useful properties of the distance func- 
tion d = dist( ■ , K) from K. We denote by Y e and K e the level sets d = e in 
P and M, respectively Thus, these level sets are equidistant respectively 
from T and K. It is immediate that K e is the Killing cylinder over T e . Since 
T is assumed to be C 2 ' Q , the function d is also C 2 ' a at points of *$?(R x Q e ), 
where Q € C is a small tubular neighborhood of T. Thus, we may define 
Fermi coordinates on $(1 x Q e ) as follows: for ge f(Rx Q e ) we associate 
coordinates s\ d by q = exp p d rj when p = p(s, s 1 , . . . , s n_1 ) in if. Then 

|W| = I. (2) 

From this it follows that 

d'di-j = 0, (3) 
where d % = g lj dj as usual. We also have 

(V dd Vd,d d ) = ^d d \Vd\ 2 = 0. 
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Therefore, 



Ad\ d=t = (V dd Vd,d d ) + f(V ds Vd,d s ) +a^(V 9i Vd,d j ) 

However Vd\d= e = d d is the unit inward normal vector field r] e to the 
equidistant cylinders K € . Denoting 

« e = /(V a .9 a ,Vd>, 
we have at points of K e that 

Ad = -K e - a'-'hijU ) = -K e - (n- l)h e , (4) 

where 6y(e) and h e are the second fundamental form and the mean cur- 
vature, respectively, of the hypersurface r e relatively to the unit inward 
normal r} e . It follows that 

where H cy \(e) is the mean curvature of K e with respect to r/ e . Its Weingarten 
operator is denoted by A e . The mean curvature of T and K are denoted, 
respectively, by h and H cy \. 

Remark 1. All of the above calculations on the distance function remain 
valid if we replace Vt t by the larger subset VLq in Vt consisting of the points 
which can be joined to Y by a unique minimizing geodesic. It was shown 
in lllOl that in this set the function d has the same regularity as T. 

Throughout this paper, the ambient Ricci tensor in a given direction v 
is defined by 



Ric M (^) = y](R(ei,v)v,ei), 



where R is the curvature tensor in M and e^, . . . , e n , v is an orthonormal 
basis. We follow [71 or ilTTI and use the fact referred to in Remark [TJ on the 
distance function in fi for proving in terms of the notation we fixed above 
the following result. 
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Lemma 1. Assume that the Ricci curvature satisfies Ric M > — ninf r H 2 V Let 
yo G r be the closest point to a given point x G T e c ri . Tfoen, we foroe 

Proof: We use local coordinates s° = s, s 1 , . . . , s™" 1 as defined above. At 
d = e and since d d is the unit speed of a geodesic, we have 

- ^(AAAj) = d d \ d=e (V 9i d d ,dj) = (V dd V dt d d) d 3 ) + (V di d d ,^aA) 

= -(Rid,, d d )d d , dj) + (AA, A e dj) (5) 
where R is the curvature tensor of M. On the other hand, 

^-(A e di, dj) = (V dd AA, dj) + (AA, ^a d dj) 

= ((V dd A e )d u dj) + (A t d^ V dd di) + (AA, V 9j dd) 

= (A'A,dj)-2(AA,AA)- (6) 

From (O and © we obtain the well-known Ricatti equation 

A' e - A\ - R e = 
where R e = (R(; d d )d d , -) | rf=e . Taking traces we obtain 

n^H cyl (e) = 3 d trA t = trV dd A e = tr(A 2 e + R e ) > nH^(e) + Ric M (d d ). 

From our hypothesis on RicAf we have that z = H cy \(d) — inf r H cy \ satisfies 

z'(d) > z 2 (d) - z 2 (0) = (z(d) + z(0))(z(d) - z(0)). 

Thus z'(d) > c(z(d) — z(0)) in some interval d G [0, d > 0] for a constant 
c > 0. It follows easily that H cy \(t) does not decrease with increasing d. | 

2.2 The mean curvature equation 

In what follows, we assume that the mean curvature H of the Killing graph 
E is a function on 0,. Computing at q = ^(u(p),p) G S, we have 

nH[p) = -tr s ViV = -div s iV. 
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Hence, if ei, . . . , e n is an orthonormal tangent frame at q in E, then 

-ntf = ^(V ei iV, 6i ) = ^(V e ,iV, e,) + (V N N, N) = div M ^-(q). 

i i 

Consider a normal coordinate frame d\ , . . . , d n with = 5{j at p e P. 
Then, we have an orthonormal frame i?o = ^^dg, E^q) = ^ <9j(p) at q. 
Using this frame and (H), the divergence in the formula above becomes 

,. V$ rat f \ f A . o ,/v7 ^* Vm «\ /v7 ^* Vm p\ 

dlVM H^ = < V (f^'^ + w M s ~ f ^ 9s ~^T 7 s ^ ~ ^ Vi?!_ ^' ^ 

The Killing equation implies that K = d s is divergence-free, and that / and 
W do not depend on s. Thus, 

j 

i 

= /(-^>V 9s 9 s )-^(V a , — ,0,). 

i 

Since ^ preserves the field <9 S and / is constant along the flow lines, we 
have 

From these calculations it results that 

i 

where the expressions on both sides are now evaluated at p e P. Since the 
di's are orthonormal at p and P is totally geodesic, we may write 

div P (^) - f(^, V 9s d s ) ~ nH = 0. (7) 

The Killing equation implies that the field Vg s d s is tangent to the leaf P. 
On the other hand, it is easy to see that 

(V/,Vn) = 2/ 2 (Va s 9 s ,Vn). 
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Using this expression and after some manipulation, one proves that an- 
other way to write out © is 

w^ 3 ~ ~ w^ if + w2){Vu > ^ dA) ~ nH = (8) 

where u^j is the Hessian of u in terms if the coordinates x % in P. Denoting 



and 



««(«,V«) = ^-^) (9) 

b(x, v«) = ~u + w 2 ){\/u, fv B ,a,), (io) 

the mean curvature equation becomes 

Q[u] = + b - nH = 0. 

The matrix a 1 ^ is positive-definite with eigenvalues 

A = and A = — 

W 3 W 

with multiplicities 1 and n — 1 corresponding to the directions parallel 
and orthogonal to Vu, respectively. Notice that A < A since / < W 2 by 
definition. 

Let be a C 2,a function on T. The Killing graph of is a codimension 
two submanifold of M. Thus, S is a Killing graph with prescribed mean 
curvature H and prescribed boundary given by the graph of <p if and only 
if u solves the Dirichlet problem 

Q[ u ] = 0, u\ r = (11) 

for a quasilinear elliptic PDE. We may apply maximum and comparison 
principles to ((TT) . Indeed, this follows from our hypothesis that the func- 
tion H does not depend on u (see, e.g., 0, Ch. 10). 



3 Height estimates 

In this section, we obtain apriori C° estimates for solutions of the Dirichlet 
problem ((TT]) . We divided the exposition in two cases concerning different 
assumptions on the ambient Ricci curvature. 

From now on, the distance function d is regarded as the distance from 
T on the totally geodesic hypersurface P. 
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3.1 Killing cylinders as barriers 

We assume that the ambient Ricci curvature satisfies 

Ric A / > -n inf H 2 yl , 

i.e., the hypothesis of Lemma 1. In this case, we construct barriers for u in 
(|TT|) of the form 

ip(x) = sup + h(d(x)) (12) 
r 

where the real function h will be chosen later. Along f2 we have 

tpi = h'di and ifi-j = h"didj + h'di-j. (13) 

As in and © above we have | Vd\ 2 = d l di = 1 and d l di.j = 0. Moreover, 
((D now reads as 

cf ;i = cr^djy = —{n - l)/i e . 
It is convenient to write W 2 = f + h' 2 . Then © yields 

Q[(p]+nH = ^-^^)-^U + W 2 )UV d A^) 

= w(~ h ' {n - l)K + h "- h ^)-w u+ w2) (/ ^' 

= ^(h" - h'(fV d3 d s ,Ve)) ~ ^((n~ l)K + (fV ds d s ,Ve))- 

We choose for ((12"]) the test function 

P CA 

where A > diam(fi) and C is a positive constant to be chosen later. Then, 

h > = e c(A-d) and h » = _ch'. 

Since the mean curvature of the equidistant cylinder K e is given by 

nif cy i(e) = (n - l)h e + k 6 , 

we get 



fh' h' 



12 



Assuming sup n \H\ < inf r H cy \ and using Lemma[T} we obtain 

fh' h' 
Q[ip}+nH<-l^(C + K e )- w n\H\. 

Observe that f /W 2 < 1. Moreover, as C — > oo we have that 

h! _ ti 

Choosing C»0 such that C + k € > 0, we obtain 

Q[y>] < -n(H + \H\) < 0. 

We conclude that at points of f2 it holds that 

Q[(p) < Q[u) = 0, 
ip\r > u|r- 

We now prove that cp > it on 0. By contradiction, assume that there 
exist points for which the continuous function u := u — if satisfies u > 0. 
Hence m := u(y) > at a maximum point y G f2 of u. Choose a minimizing 
geodesic 7 joining y to T for which the distance d = d(y, T) is attained. 
Thus, j(t) = exp yo tr), < t < d, starts from a point y G T with unit 
speed r]. Since 7 is minimizing, we have d( 7 (i), T) = t and the function 93 
restricted to 7 is differentiable with <^'( 7 (i)) = e c ' < - A_ *' ) . Since the maximum 
of u restricted to 7 occurs at t = d, i.e., at the point y, one has that 

u'( 1 (d))-< P '( 1 (d))=u>( 1 (d))>0. 

This implies that 

(Vu(y), 1 '(d)) > ^( 7 (d)) = e C(A " d) > 0. 
In particular Vu(y) 0, and hence the level hypersurface 

S 1 = {a; G ft n S r (y) : u(x) = u(y)} 
is regular for small radius r. Along S we have 

u(a;) + <p(x) = u(y) + tp(y) > u(x) + <p(y), 
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and since <p is an increasing function of d( ■ , Y) it follows that d(x, Y) > 
d(y, T) = d. From this we conclude that the points in S are at a distance at 
least d from Y. Since S is C 2 it satisfies the interior sphere condition: there 
exists a small ball B E (z) touching S at y contained in the side to which 
Vu(y) and "/{d) points. Thus, the points of B £ (z) satisfy u(x) > u(y), and 
hence 

(p(x) + m > u{x) > u{y) = (p(y) +m, x G B £ (z), 

where in the first inequality we used the definition of m. Again because cp 
is an increasing function of d, we have d(x, Y) > d on B £ (z) and therefore 
this ball is contained in the interior of f2 far away from Y. This allows us to 
extend the geodesic 7 through B £ (z). We claim that the center z of the ball 
is contained in this extension. Otherwise, the broken line consisting of 7 
and of the radius in B £ (z) from z to y has length smaller than a minimizing 
geodesic joining z to y G T (for a suitable small e such a geodesic must 
cross the level hypersurface S at a point x ^ y at distance to Y greater than 
d). Thus, if there exists at least two distinct minimizing geodesies joining y 
to T, then the point z is contained in the extension of both geodesies after 
its intersection at y. Choosing e sufficiently small, we see that this config- 
uration is not possible (the construction we made above applies to both 
geodesies). This contradiction implies that the maximum point y belongs 
to VLq. However, in this case, u(y) < 0, a contradiction. We conclude that 
u < ip throughout Vt and therefore cp is a continuous super-solution for the 



Dirichlet problem (|TT|) . 



In a similar way, we may construct lower barriers for u, that is, contin- 
uous sub-solutions for ((TT) . It is clear that the existence of these barriers 
implies the desired C° apriori estimates. 



3.2 Geodesic spheres as barriers 

Next we assume only that the Ricci curvature has a finite lower bound, 
that is, 

RicAf > ~{n - l)k, 

for some positive constant k. In this case, we present a strategy for ob- 
taining height estimates for u. Our method relies in a Hessian comparison 
theorem and takes geodesic spheres as barriers. 

We fix a point p G P and consider the function r = dist(po, • )• Let 
M n+1 (— k) be the hyperbolic space form with constant sectional curvature 
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—k and rh yp a distance function on it. The Laplacian Comparison Theorem 
(cf. Ifl5l , p. 5, Corollary 1.1) yields 



Ar < A f 



L (-fc) r hyp 



at corresponding (equidistant) points, whenever r is differentiable. This 
implies that if we consider geodesic balls B in M (outside the cut locus of 
M) and -Bhyp in EP +1 (— k) with same radius r , then the mean curvatures 
of the respective geodesic spheres calculated with respect to the gradient 
of the distances satisfy 

-H dB < -H dBhyv = -Vfccoth y/kr . 
Thus, if we assume that our function H satisfies 



1 jmaxnl-Hl 
r < — = cotn — — 

Vk Vk 



then we have that 

\H\ < Hg B . 

Thus, if we suppose that the domain Q is contained in a geodesic disc of 
radius r in P, then the corresponding geodesic sphere in M is a barrier 
for C° estimates for the problem (|TT]) . Indeed, it suffices to move such a 
sphere along the flow lines of Y and then apply the maximum principle at 
a first tangency point between the graph S and the moving spheres. 

Remark 2. For constant mean curvature H it was shown [6] that an apriori 
height estimate exists if Ricm > —nH 2 . This is achieved by constructing a 
function that is subharmonic but only if H is constant. Up to this estimate, 
all the aforementioned results on Killing graphs in the Introduction for 
different ambient spaces follow from Theorems [1] and |2] in this paper. 



3.3 CMC spheres as barriers 

Next we assume that the induced metric in P is rotationally invariant. 
More precisely, the metric in P is of the form 

dr 2 + £ 2 (r)d6 2 
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in terms of coordinates (r, 6) G R + x §™ _1 , where d9 2 denotes the usual 
metric in §™ _1 . We also assume that g = g(r), that is, the norm of the 
Killing field does not depend on 9. In this case, the ambient metric is 
written in terms of cylindrical coordinates s, r, 9 as 

g 2 (r)ds 2 + dr 2 + £ 2 (r)d9 2 . 

and M is a doubly-warped product with respect to warping functions of 
the coordinate r. 

A rotationally invariant hypersurf ace S is parametrized by an immer- 
sion R x §™ _1 — > M whose coordinate expression is 

(u,9)^(s(u),r(u),9), 

where u is the arc-length parameter of the profile curve 9 = 9 . This means 
that u is defined by g 2 s 2 + r 2 = 1 and that the induced metric in S is 

du 2 + i 2 {u)A9 2 . 

If S has constant mean curvature H , then it satisfies a first order equa- 
tion given by the flux formula. In the flux formula in the Appendix, we 
put T as the geodesic circle of radius r — r{u) which is the intersection of 
S and the leaf F s , where s = s(u). Thus we consider D as the geodesic 
disc in ¥ s with radius r = r{u). The co-normal v is the unit velocity vector 
sd s + rd r . Finally, the Killing vector field Y corresponds to the coordinate 
vector field d s . Thus, one has 

(Y, v) = g 2 (r(u))s and (Y, N D ) = (Y, = g(r(u)). 
Hence, plugging these expressions in the flux formula gives 

2 it I I „t»-lj„Jfl i / -„2/-n-l 



c = nH g+ sg z = nH / gC~ drd9 + / sg 2 C~ d6. 
Jd Jy Jo Js"- 1 is"- 1 

Since we are integrating at a fixed value of s and therefore at a fixed value 
of u, we have 

nH [ gC^dr + s(r( M ))f ? 2 (r)e"- 1 (r) = — . 

JO 
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The derivative of this expression with respect to u gives a second order 
ODE which characterizes CMC rotationally invariant hypersurfaces in M. 
Compact solutions satisfy r = at the points of maximum and minimum 
height, where s = 0. Hence, these compact examples correspond to take 

c = 0. 

On the other hand, at maximum points for r we have g 2 s 2 = 1 and 
therefore 

pro 

nH / g(r)C-\r)dr + ^ro^fo) = 0, 



where r is the maximum value of r. Thus the mean curvature of compact 
rotational examples with maximum radius r satisfy 

Jo £?( r )£ [r)dr 

If we assume that O c P is contained in a geodesic disc with radius r , 
then these CMC spheres are barriers for the height of a Killing graph with 
prescribed mean curvature H(x) satisfying 

n\H(x)\ < F(r ). 

The value of F(r ) may be explicitly given in particular cases such as 
M n+1 (—k) and W l (—k) x R. 



4 Boundary gradient estimates 

Our task now is to produce apriori gradient estimates for the Dirichlet 
problem ((TT) . In order to do that, we use barriers of the form w + along a 
tubular neighborhood fi e of T as defined in Section 2.1. Here, w = ip(d(x)) 
for some real function ip to be chosen later and d = dist( ■ , T). Moreover, 
the boundary data <fi was extended to f2 e by 4>(s\ d) = 4>(s l ) for simplicity. 
A simple estimate gives 

Q[ w + = a ij (x, Vw + V<p)(w i;j + (f) Mj ) + b(x, Vw + V0) - nH 
< a ij Wi.j + A|0| 2 , Q + b - nH, 

where a* J and b are given by © and (|T0l) . Thus, 
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On one hand, we deduce from (O and © the expressions 

wViUiy = (^') 2 d i d j {^"d i d j +^'d i . j ) = (^) 2 ^"|Vrf| 4 = (V>')V 

and 

^ywy = ip'dtftWdidj + ip'di.j) = VV'I Vd| 2 *>j = W(Vd, V0) = 
and 

o'o '//•,., = <f) i (j) j (ilj"d i d j +ip'd i , j ) = ^"(Vd, V0) 2 -0W&ij(e) =-^'<f> i <ff i b ij {e 
In particular, we obtain 

^-^w«* = £(i-^!) = £(/ + |v#i») 

since W 2 = / + (^') 2 + |V0| 2 and 

On the other hand, we deduce from © that 

Aw = V" + ^'Arf = V" - (n - 1)^6- 

We conclude that 

= -^(" - * + ^(/ + IV0I 2 ) + ^V*M<0- 
A suitable expression for b is 
1 



W 3 



U + W 2 ){Vw + V<PJV da d s ) 



since Vw = and k £ = f(Vd s d s ,r) e ). Therefore, 

Q[w + 0] < -|((n - l)K + « e ) - n# + A|0| 2 , Q - ^ 

-tMw* + l)</Vo.a., V0) + £(/ + |V0| 2 ) + ^V6y(e). 
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Finally, using that A = 1/W, we get 
W 3 Q[w + <j)}< -mP'H cyl (e)W 2 - nHW 3 + \(j)\ 2 , a W 2 - ip'f Ke 

-(fVaA, V0)^ 2 + VWM<0 - f 2 (v 9s d s , v^) + ru + 1 V0| 2 ). 

Now define 

ip(d) = fi\n(l + Kd) 
for certain positive constants fj, and K to be chosen later. We have 

^' = T^7i and r = --m 2 - 

1 + Kd jj 
We choose /i in such a way that — ► as K — > oo. It suffices to take 

C 

/i 



ln(l + X) 



for some positive constant C to be chosen later. In this case, as K — ► oo 
one has 

It also holds that ^ ~ 1 as if — > oo. Thus, at points of T the last inequality 
(asymptotically) becomes 

iy 3 Q[^ + 0] < -n^i + W 3 - -(/ + |V0| 2 )(^') 2 

+ (|0| 2 ,a - (/V 9s 9 s , V0»(V/) 2 - (fK + WbiM - f 2 (V ds d s ,V<P) 

Therefore, assuming that H cy \ + H > and choosing if large enough, we 
assure that Q[w + 4>] < on a small tubular neighborhood fi e of T and that 
w+4> > w|r e +0onboth boundary components. Therefore, w+4> is a locally 
defined upper barrier for the Dirichlet problem (|TTT) . A lower barrier may 
be constructed in a similar way. 



5 Interior gradient estimates 

The last step in providing apriori estimates for (ITT) is to verify that Vu 
satisfies a kind of maximum principle for a third order equation obtained 
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from differentiating Q[u] =0 and contracting the resulting equation with 
the gradient itself. 

Using a suitable test function taken from [19] and Ricci identities allows 
us to eliminate third derivatives and to obtain global estimates for \Vu\ in 
terms of the height and boundary C 1 estimates. 

We suppose momentarily that u E C 3 (f2). The usual regularity theo- 
rems guarantee that the estimates we will obtain are also true for a C 2,a 
function (see [7j and lfi~8l ). Equation © may be written as 

■v »r.j r>. (14) 

where 

A ij (x, Vu) = W 2 a ij - uV (15) 

and 

B(x, Vm) = (/ + W 2 ) (fV ds d s , Vu) + nHW 3 . (16) 
Differentiating covariantly with respect to the metric cr^ on P yields 

V k W 2 = V fe / + Vk{w>Uj) = fk + u 3 . k Uj + u 3 Uj- k = f k + 2u J Uj- k . 



We obtain from CL4J) that 

W k B = ((/* + 2u l Ul .k)cr lJ + W 2 W k a 11 - Vi(uV'))«y + A ij u i;jk . 
Contracting with u k gives 

u k V k B = (f k u k + 2u j u k u j . k )u i . i - -In' ti, :j ti{ k ti k ' + A ij u k Ui- jk . (17) 



Following [19J we define the function 



r = e 2Cu v, 



where v = u l Ui is the squared norm of the gradient of u and C > a 
constant to be chosen later. Then, 

n = e 2Cu {2CvUi + Vi). 

However, 

= Vi{u 3 Uj) = u 3 .^ + vPuj ; i = 2u^Uj-i. 
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Thus, 



n = 2e 2Cu (Cv Ui + u j Uj ,i). (18) 



If r achieves its maximum on Y then we have a bound for v in Vt as desired. 
Hence, we may assume that the maximum is attained at an interior point 
Pq EVL where we have 

v?v,j. t i = —Cvui. (19) 

Differentiating ((IB"]) yields 

T i .j = 2e 2Cu (2C 2 u i u j v+2Cu j u k u k - i i+Cui- d v+2Cu i u k (20) 

From the maximum point criterion we have r i;i < 0, and from the elliptic- 
ity of CEU that 

A ij r Mj < 0. (21) 
From ((T9l> we obtain atp that UjU 3 ^ = u j uj-i = —Cuiv. Therefore, 

//•'// ;: ,//' = -Cv 2 (22) 

and 

u?u 3 .^\ h u k = C 2 v 2 u k u k = CV. (23) 

We may assume that Vu is non-singular in a neighborhood of p . Other- 
wise, we are done. Hence, we choose local coordinates at p asking x 1 to 
parametrize the trajectories of Vu and such that 

di = Vu/\Vu\ 

along the trajectory at p . The remaining coordinates parametrize the level 
sets of u and are assumed to be orthonormal and geodesic at p . Thus, we 
have at p that 



Ui 



\Vu\= v 1 ' 2 and u, = for j ^ 1. (24) 



We also have = u\ = Uji. Then (|T9|) gives 

iii-i = u n = —Cv and = -u^ = for z ^ 1. (25) 

If necessary, we rotate the coordinates to assure that Uij is diagonal at p . 
Then (fl4|> at p becomes 

(/ + v)Au = (/ + = 5 - mt u = 5 - Cv 2 . (26) 
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We now use the Ricci identities for the Hessian of u 

^i;jk ^k;ij Rijkm^ i 

where R is the curvature tensor in P. Thus, we obtain that 

A^^u^k = ((/ + v)a ij - u i u j )u k u i;jk 

= ((/ + v)av - uV> fe (u fc;ii + i^ fcm « m ) (27) 
= .V^a'a,,,, + {f + v)a^R ijkm u k u m 
since Rij km u l u^u k u m = 0. Moreover, we have from ([22]) and ((261) that 

(/ fc u* + 2u%S ;fc )^ = (/ fc u fc - 2Cv 2 )D 1 (28) 
where L>i =: (B - Cv 2 )/(f + v). From CO, d23]>, © and (28) we deduce 



u fc V fcJ B = (/ fcU fc -2Ci; 2 )D 1 -2CV+^u fc u fc . ii + (/+i;)^'^u fc w m . (29) 
We conclude that 

.l'V'»/„;. ; = -f k u k D 1 -(f^)a i ^R ij ^u k u m +2C^v 2 f+2Cv 2 D 2 -^ k V k B (30) 

where we denote D 2 = (B — Cf 2 )/(f + v). 

Next we write down the maximality condition (|2T) for r at the point p . 
A straightforward computation using (fl5]> , (|24]) and (|25l> in (|20"1) yields 

> A ij (2C 2 UiUjV + 2CujU k u k -i + Cui-jV + 2CuiU k u k -j + u k jU k;i + u k u k;ij ) 
= (/ + v)(Cvuu + u^Utf) - 2CV/ + ^V/V.y, 

Using that Ujjiijj > C 2 v 2 and (|26|) we obtain 

A^V/Vjy < CV/ - CBv. (31) 

Combining the expressions (|30|) and (|3l) we get 

-f 1 v 1 ' 2 D 1 + (/ + v)a ii R ik j m u k u m + CV/ + 2CV 2 D 2 + C5w + w 1 / 2 ViB < 0. 

Let i? be the minimum eigenvalue of the Ricci tensor in P in the direc- 
tion of Vu. Thus, we have Rv < a lJ R ik j m u k u m and multiplying the above 
expression by / + v yields 

Rv{f+v) 2 +C 2 v 2 f{v-f)+CBv{f+3v)+f 1 v 1/2 {Cv 2 -B) + {f+v)v 1 / 2 V 1 B < 0. 

(32) 
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Next we compute the last term on the left hand side of (|32|) . From (|T6|) we 
obtain 

5(2, V«, v) = (2f 2 + fv)T 00 ^ + nH(f + v) 3/2 (33) 
where T 00j = (JijT l m . Differentiating B = B{x, Vu, v) we have 

V k B = B k + B u] u\ k + B v v k . 

We obtain at p that 

u^ViB = Sit; 1 / 2 - CB u ,v 3 ' 2 - 2B v Cv 2 . 

Replacing this in (|32~1> yields 

+ v) 2 + CV/(t; -/) + (/ + v)B l v 1 ' 2 + hv lj \Cv 2 - B) 
+Cv((f + 3v)B - (f + v)(B u iv^ 2 + 2B v v)) < 0. 



Using (|33|) we have at p that 



Bx = (4f + v)f 1 Tv 1 / 2 + (2f 2 + fv)V 1 rv 1 / 2 + nH 1 (f + v) 3 / 2 + ^nH(f + v) 1 / 2 f 1 
and 

S tt i = (2/ 2 + fv)T 

and 

s v = /r^ + 3 ni?(/ + u) i/2 j 

where T =: r 001 . In particular, 

5 1 t; 1 / 2 = 2/(2/ 1 r+/V 1 r)^+(/ 1 r+/V 1 r)i; 2 +n(/+t;) 1 / 2 t; 1 / 2 (iJ 1 (/ + ,;) + ^ /l ) 
and 

(/ + 3v)B - (/ + v)(B u .v 1 ' 2 + 2B v v) = 2f 2 Yv 3 ' 2 + nHf(f + vf' 2 . 

We may assume that max^ f < v since, otherwise, we are done. Then, we 
have f + v <2v and 

B lV 1 ' 2 > {AfhT + 2f 2 V l T-(Z/V2)n\Hf l \)v + (f l T + fV x T - 2 3 / 2 n\H 1 \)v 2 
= A x v + A 2 v 2 . 
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Moreover, we have 

(/ + 3v)B - (/ + *;)(iW /2 + 2B v v) > (2f 2 T - 2 3 / 2 n\H\f)v 3 ' 2 = A^v 3 ' 2 
and 

hv ll2 {Cv 2 -B) = f 1 v 1/2 (Cv 2 -(2f 2 + fv)Tv 1 / 2 -nH(f + v) 3/2 ) 

> -h{2f + fv)Tv - 2 3 / 2 n\Hf 1 \v 2 + Cf^ 2 

= Civ + C 2 v 2 + C 'hv b / 2 . 

Replacing in ([34]) and grouping equal powers of v , we obtain 

(C 2 f + R + A 2 )v 3 + C(A 3 + f x )vV 2 + (2Rf - C 2 f 2 + A, + fA 2 + C 2 )v 2 
+(Rf 2 + A 1 f + C 1 )v<0. 

Notice that all constants in the expression above depend only on H and Vt. 
We choose C large enough so that the coefficient of v 3 in the left hand side 
is positive. With this choice we conclude that v < C atp for some constant 
C depending on the height estimates for u and on the data H and Vl. 

6 Proof of the theorems 

We apply the well-known continuity method to the family of Dirichlet 
problems 

Q a [u] = 0, u\ T = cr0, 

where a G [0, 1] and 

Q a [u] = a^Ui-j + b — naH. 

The subset / of [0, 1] consisting of values of a for which the above Dirichlet 
problem has a C 2,a solution is non-empty, since G I. The hypothesis 
H u = implies that / is open. This follows from a standard application of 
the implicit function theorem. The closedness of / follows from the apriori 
estimates we had proved. Thus, the continuity method assures that 1 G I. 

In order to prove uniqueness, it suffices to reproduce the proof pre- 
sented in |6l. ■ 

We point out that our existence results still hold if 4> is only assumed 
continuous. We may approximate 4> uniformly by smooth boundary data 
and use the interior gradient estimate to obtain strong convergence on 
compact subsets of Vt. A local barrier argument shows that the limiting 
solutions achieves the given boundary data. 
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7 Appendix: a flux formula 



Consider a hypersurface £ in M with boundary T and another hypersur- 
face D such that MUD bounds a domain U in M. Let be a function 
that is constant along the flow lines of Y. Choose an orientation onMUD 
given by unit normal vector fields N along £ and N D along D. By the 
Killing equation, we have 

div M HY = Hdiv M Y + (VH, Y) = 0. 

Thus, applying Stokes theorem to the domain U, one has 

0= f div M HY= [ H(Y,N}+ ! H(Y, N D ). 
Ju Jt, Jd 

However, if we consider an orthonormal adapted frame N, e±, . . . , e n to E 
we have from the Killing equation 

o = ]T(v ei r, ei ) = $](v ei r T , + ^(v 6l iv, ei ) = div s y T - nH{Y, n). 

i i i 

Thus, applying Stokes theorem to £ we obtain 

njH{Y,N) = J{Y,u) 

where v is the exterior unit co-normal of E along T. The two expressions 
we obtained above may be gathered at the flux formula 

n [ H(Y,N D }+ [ (Y,u}=0. 
Jd Jr 

An useful variant of this reasoning is to consider a hypersurface T in E not 
homologous to zero in E such that T = dD. Then we have that 

njH{Y,N D )+ J{Y,u) = c, 

where c is a constant depending only on the homology class of T. 
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